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ABSTRACT 

PLATE ANALOGY OF BEAK GRILLAGE FOR 
DYNAMIC ANALYSIS 

by 

rajtshwar dayal sharma 

M.Tech. (Mechanical) 

Indian Institute of Technology, Kanpur 

*-» March, 1973 

In the present work dynamic analysis of beam 
grillage has been undertaken*- The analysis is based «n 
the rational proposition of orthotropic plate analogy 
formulated on the basis of energy equivalence, replac- 
ing the Timoshenko's plate analogy which is an intutive 
proposition* 

The analogy is validated for the case of 
orthotropic grid with simply supported edges by com- 
paring the results of the frequencies with those obtained 
by Finite Element method. An attempt has been made to 
extend the proposition to the case of grillage with 
clamped edges but is limited to fundamental frequency 
«>f the orthotropic grillage. 



CHAPTER. I 


INTRODUCTION 

The use of computer for dynamic analysis of team 
grillage is very common now a days- This is time consuming 
and costly. In the present work, an attempt has been made 
to establish plate analogy of beam grillage such that the 
analysis is simplified. 

A grillage is a structure consisting of. two sets 
of orthogonal or non orthogonal beams rigidly connected at 
the intersections. Some of its applications are in building 
floors (concrete slabs embedded in grillage), ship decks, 
bridge floors, missile ground facilities, electronically 
steerable radar systems etc. Several studies have been 
made of static and dynamic analysis of beam grillage. 

These are classified into two groups. 

(a) Analogy methods 

(b) Methods using discrete approach 

Discrete approach is very suitable for analysis 

of such systems but is laborious and time consuming. 
Furthermore if cne uses consistent mass matrix approach, 
it becomes difficult to solve higher order grids due to 
limited memory space of the computer. Terms like isotropic 
and orthotropic grillages are commonly used in this disser- 
tation and are explained as follows; 
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Isotropic beam grillage : 

A grillage made of two sets of orthogonal beams, 

equally and uniformly spaced in both directions and having 

* 

same geometric and material properties, is called an iso- 
tropic beam grillage. 

Orthotropic beam grillage : 

An orthotropic grillage is made of two sets of 
orthogonal beams, each having different geometric or 
material properties, -dth uniform spacing of the beams in 
each direction but the spacing in one direction may be 
different from that in the other. 

The above terms have been defined in a manner that 
while proposing the analogy, these have similar meanings 
as those used for plates. 

1 .1 STATEMENT OF THE PROBLEM : 

The present 'fork deals "ith formulating an analogy 
for the dynamic analysis of the beam grillage for finding 
its material frequencies. As mentioned previously discrete 
analysis is time consuming and as such an attempt is to 
find an analogy of the physical system (grid). Thus the 
problem under consideration is concerned with finding a 
physical continuous system equivalent to that of beam 
grillage in m, n mode. This implies finding unknown 
parameters, the flexural rigidity and the thickness, of 
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the equivalent plate in terms of the 'mown parameters of 
the grillage. 

The statement of the problem being undertaken in 
this •■ r ork is thus 

"To establish the plate analogy of the beam 

grillage for its dynamic analysis". 

1 .2 PROPOSED ANALOGY ; 

The analysis developed is based on formulation 
of energy equivalence. Thus for dynamic analysis of a gri- 
l lag e'^; it' is said to be replaced by an equivalent conti- 
nuous system, plate, if the maximum strain and kinetic 
energies of the two systems are equal. In the present 
work, transverse vibration of the grillage is being consi- 
dered, as such, equating the maximum strain energy of the 
grillage to that of the plate in bending and equating 
maximum kinetic energy of the grillage to that of the 
plate, yield the required unknown parameters of the equi- 
valent plate in terms of tbe known parameters of the 
grillage. 

The parameters thus found, are substituted in the 
frequency equation of the equivalent plate which directly 
gives the frequencies of the given beam grillage. 

In the present work, beam grillages with the follow- 
ing two sets of boundary conditions are considered. 



1. Grillage, simply supported co the "boundaries and 

2. Grillage clamped on the boundaries. 

However the proposition is general and can be 
applied to grillage having any set of boundary conditions 
if the corresponding equivalent plate solution is known. 

In the second chapter the review of the earlier 
work is given. Third chapter deals mainly with the formu- 
lation and method of solution of the problem. Case studies 
are presented in chapter four and the results obtained from 
discrete analysis and analogy method are compared. These 
results are also compared with the results of earlier work 
done. Chapter V concludes with the validity of the analogy 


proposed. 



CHAPTER II 


LITERATURE REVIEW 

It has already been mentioned that the present 
work deals with the dynamic analysis of grillage. Several 
studies have been made on the subject and these can be 
divided in two groups; 

1 . Studies using discrete approach and 

2. Studies which replace the grillage by an 
equivalent continuum,, plate. 

In this chapter, a review of the earlier work on. 
grillages is given. 

2.1 DISCRETE APPROACH ; 

In. these methods, the system under consideration^ 
is divided into small elements such that the overall beha- 
vior of the system is unchanged. The stifftiess and mass 
properties of these elements are- -determined which aftQf 
being assembled for the whole system give its mass and 
stiffness properties. 

The earlier work on the vibration of grillages is 
known to be that of Cox and Denke^ who formulated the pro- 
blem into matrix form, suitable for use with a digital 
c omput er . 
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J.P. Ellington & H. McCallion have taken into 
account the discrete nature of the beams and considered the 
grillage system to be made up of two orthogonal sets, of 
beams with equal spacing attached at the joints so that nt 
moments, bending or torsional, are transmitted from one set 
to another. The mass of the beams is taken as lumped at 
the joints (thereby reducing the problem to finite degree 
of freedom) so that there is an equal point load at each 
internal joint. Taking the bending moments and deflections 
under each load as unknowns, the equations of continuity 
and equilibrium, taken with appropriate boundary conditions, 
furnish the required frequency equation. They used the 
finite difference technique, coupled with three moment 
equations, for developing frequency equation. 

A more general solution was presented by Thein 
Wah^3, who extended the method of J.P. Ellington and H. 
McCallion by taking distributed mass of the beams instead 
of lumping the masses at the joints. He further extended 
the method so as to investigate the effect of St. Venent 
torsional inertia on the transverse vibrations. However 
due to complexity of the formulation, the solutions obtain- 
ed are for bending deformations only. 

It is to be remembered that the solution of 
£ 

Ellington and McCallion 0 is not an approximate solution 
but rather an exact solution of an approximate formulation . 
Possibly as a consequence of this, the lowest order 
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frequencies are sometimes under estimated and sometimes 
over estimated. 

23 

The in Wah r s work was extended by Chang and 

3 

Michelson and the complicated equations were simplified 

q 

with the application of Raleigh Ritz method. Chang's 
method is specially useful for the vibration analysis of 
the ship structures. 

1 6 

Miller and Kamm employed the new technique of 
matrix formulation for the dynamic response of grid struc- 
tures by using the lumped parameter method for which the 
static spring is composed of flexural and torsional stiff- 
ness coefficients formulated from the classical slope 
deflection equation. 

If 

Cheng analysed the torsional-flexural vibrations 
by using the general dynamic matrix of distributed mass 
system for the orthogonal or non orthogonal grillages com- 
posed of either prismatic bars ‘or thl$. walled open section 
members. The differential equations are derived and the 
dynamic stiffness coefficients are th&n formulated for 
which a general computer programm is developed that is 
applicable. %o various grid works with different boundary 
conditions of fixed, free, or hinged support. The effect 
of rotary inertia, shear and bending deformation and that 
of St. Venent torsional inertia on the natural frequencies 
are taken into, account for the grid work of prismatic bars. 
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•For the grillage ^ith. wide flanged members, the additional 
effect oT warping and the associated effect of longitudinal 
and transverse -shear deformation of the flange on the vib- 
ration are also included. The effect of the sectional 

dimensions on frequencies is also considered in the method 
. k 

of Cheng which consequently leads to accurate solutions 
of higher modes. 

H- 

It should be remembered that the method of Cheng 
presents difficulties when the order of the grid becomes 
high due to limited computer memory space. 

2.2 ST J'DIES WHICH REPLACE THE GRILLAGE 
BY AN EQUIVALENT PLATE; 

For beam grillages, the idea of replacing a dis- 
crete system with an equivalent continuum has received an 
increasing attention. In these studies, the actual beam 
grid is replaced by an equivalent plate which in turn gives 
the solution of the grillage. 

Several studies have been made in which the grid 

is replaced by an equivalent orthotropic plate. For static 

2k 

analysis Timoshenko considered the grid, shown in Figure 
2.1 and assumed that if the distances a„ and b. between 
the beams are small in comparison with the dimensions a and 
b of the grid and if the flexural rigidity of each of the 
beams parallel to the x axis is equal to and that of 

the beams parallel to y axis is equal to , the plate 
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parameters are written as 
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D x =B 1 /b 1 D y= B 2 / a 1 (2* 1 ) 

He expressed the quantity D in terms of the torsional 
rigidities c^ and c^ of the beams parallel to x and y axes 
respectively as follows: 




■Oy 


M. 


3 2 


V 


y 3 - a 1 ftx5y 


( 2 . 2 ) 


Since the beam is considered as a one dimensional 
continuum, poisscn 's ratio$ is zero and the quantity 
for the grillage, defined below, 

Ex V y h 3 

D 1 " 12 (1 - V x Xly) 
is also zero. 


Thus the parameters of the equivalent orthotropic 
plate are known in terms of the known parameters of the 
given grid and the static analysis of the grid can be dene 
easily. 

22 

Using Timoshenko's results, TaraporewaHa found 
static deflection of the grillages by considering the plate 

pi* 

analogy. The relations given by Timoshenko were purely 
intutive and gave no indication of the accuracy of the 
method. Further he assumed that the spacing of the beams 
in the grid is very small, which may not be. valid in most 
of the grillages. 
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' 1 9 

Renton J.D. gave the theoretical basis for the 
relations given by Timoshenko. He took advantage of the 
resulting regularity of the grids and the fact that a large 
number of beams are used in their formation. He used the 
combination of n finite difference”, and differential calcu- 
lus, the two being related by Taylor's expansion. This 
technique gives differential equation for the deflection 
of an equivalent continuum. 


In case of planer grids, the resulting equations 
are biharmonic and are analogous to those of plates, so 
that known solutions of the plate problems can be used to 
solve grid work systems. 


For the rectangular mesh grid, shown in Figure 2.2, he 
obtained the following equation 


( A GJ-ys ^w _ 

*1 bl a 1 ; VV a-, B/ " ^ 


( 2 . 3 ) 


in which EL and GJ are the flexural and torsional 
xx x 

rigidities of beams of span in the x-direction and Eyl^_ 
and GJy are flexural and torsional rigidities of beams of 
span b-j in y-direction, q X y. is load per unit area. In all 
the cases, as the mesh sizes a^ , b^tend to zero, solutions 
tend to the corresponding plate solutions. 


Renton"' 9 m ade the same assumption as that of 
Timoshenko that is the spacing of the beams in the grid is 
small in comparision to overall dimension of the grid p-n* 



while deriving equation (2.3) he ignored terms proportion- 
al to the square of the ratio of the mesh size to the over- 
all dimensions of the grid. 

Of) 

In his further work, Benton gave 1st. and 2nd. 
order continuum approximation for various forms of epa.ce 
grids. The introduction of 2nd. order terms not only gives 
more accurate results but also indicates the range of app- 
licability of the concept of analogous continua. The 
method is inappropriate for grids when the fineness param- 
eter c = a/a-j, where a is the side of boundary and a^ is 
the length- of unit step, is less than four because the 
convergence is not assured. 

William B. Flower and Lewis C. Schmidt" analysed 
a particular form of double layer .space truss by rcpl?cm^ 
it with an equivalent plate. From consideration of the 
Partial differential equation derived from the equi librium 
equations of the truss by means of Taylor's expansion, it 
is possible to express the elastic constants for the equi- 
valent plate in terms of the grid parameters. This is 
possible if the derivatives of the deflections are bounded 
which implies that the Taylor series expansion is valid. 

Ll 

In the paper of Cheng , the results for the 

frequency parameters for a 3 x 3 square uniform orthogonal 

grid are given, which have been obtained from orthotropic 

7 

plate solution using plate analogy. The results compare 
well with the result obtained from various discrete methods. 


CHAPTER III 


mathematical formulation and method of solution 

As mentioned, the analogy is based on the propo- 
sition that plate with boundaries same as that of grid is 
dynamically equivalent to grid if its maximum kinetic and 
potential energies are same as that of grid. This equiva- 
lence gives two required unknown parameters, flexural 
rigidity, D and the thickness h of the equivalent plate, 
in terms of known parameters, elastic constant E, moment 
of the area I and given dimensions of the grillage. Once 
the parameters of the equivalent plate are known, well 
established solution of the plate can be used to fine 1 the 
natural frequencies of the grid. 

3.1 MATHEMATICAL FORMULATION : 

In the present work, analogy for the class of 
grillages is considered for which the corresponding equi- 
valent plate solutions are known. Since exact solution 
is known for the plate with simply supported edges, the 
grillages with simply supported ends are first considered. 
Furthermore approximate solutions of the plate clamped m 
the edges is also known, as such having established the 
equivalence for grillage with simply supported edges, the 
analysis has been extended to grillages clamped on the edges. 



14 


3*1*1 Assumptions : 

The grillag-e~ia._analy.sed for transverse vibra- 
tions wherein any element Is considered to be. in bending. 

»v 

The analysis is based on certain assumptions, which are 
summarized below: 

1 . The grillage is torsionless. 

2. The rotatory inertia effects and the shear effect# in 
the transverse vibrations of the grillage are 
neglected. 

3. The transverse deflections of the grillage are 
considered to be small. 

4. The equivalent plate follows Kirch off's hypothesis 
which implies'. 

(a) A normal which is perpendicular to the middle 
plane remains normal to the def ormed micldle 
surface and does not change its length. 

(b) Middle plane of the plate is riot strained. 

(c) Normal stress transverse to the middle plane ean 
be neglected. 

The analogy proposed implies that an orthotropic 
rectangular grid is replaced by an equivalent orthotropic 
plate as such an isotropic rectangular grid is replaced 



by an equivalent isotropic plate. In particular, a sqc~r- 
isotropic grid is to be considered by its equivalent i;, 
tropic square plate. 

Consider first the grillage with simply supported 

ends. 


3*1.2 Grillage simply supported on the boundaries : 

Consider an orthotropic beam grid work shown in 
figure 3*1 a, of length a and breadth b with number of beams 
N x and N y and with flexural rigidities E X I X and EyI , in 
x and y direction respectively. According to the propose'' 5 
plate analogy, the grid can be replaced by an equivalent 
orthotropic plate with same set of boundary conditions, 
refer figure 3*1 -b . To obtain energy expressions for the 
.Vibrating elastic system, its model functions are refynfed 
'To start with, it is assumed that the model functions of 
< the grid are same as those of the equivalent plate. It is 
well known that (m,n)tb mode eigen -functions of the plate, 
simply supported on the edges, is given by 

w(x,y,t) = B m Sin SZH S in Sin P mn t (3.1) 

where B mn is the amplitude and 

P^j! is natural frequency in m, n mode. 

It can be established that mode shapes of the grid found 
by finite element method fit in closely with ^odal func- 
tions of the plate at all the discrete points. As suck 








17 


the mode shape of a beam of the grid can be. adopted as that 
obtained from the modal function of the plate correspondir : 
to its discrete value of x or y. One can thus write the 
expression for the mode shape of ith beam lying along 
x -direction as 


w. 


± ( x » 


bi 

N x+ 1 


t) = 


Sm 

mn 


a 


iN x+ ' 


(3.2a) 


where , i = 1 ,2, ,N_ 

or 

tf i ^ = B « Sln Up Sin S~! Sln Pmnt 0.2b) 

Similarly, the mode shape for the ;jth beam lying along 
y-direction can be written as: 


w j = 3m Sin Sln t? S ta v* (3<3) 

X 

where j = 1 , 2, . ., N y 

i 

The grillages, simply supported on the edges, can 
be classified as orthotropic and isotropic and are discu- 
ssed separately. 


A- Orthotropic grid : 


An orthotropic grid, as defined earlier, is made 
of two sets of orthogonal beams, each set differing from 
the other in its material or geometric properties. 

Whereas an orthotropic material has two preferred direc- 
tions, the material properties remain the same if the 
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directions are reversed. For such materials, the linear 
stress -strain relation is written &s 

<s = fs} e (3.^) 


where <r is the stress vector 
6 is the strain vector 

and K) is a square matrix having 9 components, 
known as elastic constants. For plane stress problem 
matrix PU reduces to : 




Ei 2 


0 


E 1 2 
e 22 
0 


ol 


Note; en gi neer i n S practice, material cn;.„ 

tants used are E x , E^., 5^. and G. These constants are 

related to the elements of fs] in the following form ; 


E. 


E, 


'11 


1 - \)y 


)E 22 


& 


, E 1 2 = E 21 - E 22 V x - E 1 . V 


1 - v 7 

G 1 = G . (3-6) 

The governing equation for free transverse vibra- 
tions of an orthotropic plate is 


D Z-F + 2H -KX, + D + t h 2% = 0 


(3.7) 


where H = + 20.^ and 


E^h 3 

D _ = -LI— , D > = 

12 7 12 


E 22 h3 


x 


E 1p h 3 E 21 h3 

t Vj \ : - j 


12 


12 


D =2^ 
xj 12 


(3.8) 
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Let D q , the flexural rigidity parameter for the 
orthotropic plate, S.S. on the edges, be defined as 


d. = £. d + 

0 s? x s>2h2 


H + 


n 


h 


D. 


a^b' 


^ U 7 


( 3 * 9 ) 


O 

then the frequency equation for the plate is written as; 



where is the weight density and h is the thickness of 
the plate. It is necessary that for finding the frequencies 
of the given orthotropic grillage, the flexural rigidity 
parameter D Q , and the thickness h of the equivalent ortho- 
tropic plate should be known. To do so proposed plate ana- 
logy is used. 

( i ) Energy expressions for the orthotropic 
grid work : 


The potential and kinetic energy expressions Ug 
and Kg for the grid work vibrating in m,n mode derived in 
Appendix -A, are given as : 


1 


if 


U — t? B -rr Sin p "b 
g 8 mn 71 mn 


^ Ms (N x+ 1 ) 



( 3 . 11 ) 



p t 
mn 


[v (H x 


+D 


m y b(N- 



( 3 - 12 ) 


and 
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(ii) 


for the orthotropic plate; 


Similarly the energy expressions Up and Xp for 
the equivalent orthotropic plate in m,n mode are given 


below, refer Appendix -A, 


1 = l 4 p B in TT 4 Sin2 p^t 




m^n^ 


(3.1 3. a) 


which in terms of D , defined earlier, is 


1 = l B i A P D o * 4 si" 2 p»‘ 


(3.13-b) 


K p - g B mn \ M p Cos2 Pmn 1 " 


P P 


(3*14-) 


According to the proposed analogy, the maximum 
strain and kinetic energy expressions of the plate are to 
be equated with the maximum strain and kinetic energy 
expressions of the given grillage, thus yielding the 
required unknown parameters of the orthotropic plate. Thus 
equating U , eq. 3.11, with T J p , eq. 3.13, one gets 

D ° = < v- 1 : ^ ^ ( v 1 )nlf 3 

Similarly equating the kinetic energy of the orthotropic 
grid Kg, eq. 3*1 2, to that of the equivalent orthotropic 
plate Kp eq. 3 *14- , one obtains the expression 


MpAp = m x a (N^+1 ) + ra y b (N +1 ) 


(3.16) 



Since mass per unit area of the plate is 

M P = JV 

h is obtained from 3*16 as 

h = [> a (Nx+1 5 + b (I!y+1 \ 
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(3.17) 


Thus the parameters of the equivalent orthotropic 
plate D q and h, are known in terms of the grid parameters 
m x’ m y’ N y’ ^ x ’ a?b,E x and which in turn can be direct- 
ly used to get the frequencies of the orthotropic grid, 
eq. 3*10. The following simplified cases of the ortho- 
tropic grillage are considered : 

(a) Orthotropic grid having beams equally spaced in 
both directions. 


(b) Square orthotropic grid having beams equally spaced 
in both directions. 


( a ) Qrthotropic grid having same spacing in both 
directions : 


Having found the flexural rigidity parameters D Q 
and thickness h of the equivalent orthotropic plate, it is 
easy to get these for special cases. For 


a _ b 

Ny+1 Fjc : i-l 


(3.18) 


Using equations (3*15), (3*17) and (3*18) one obtains the 
two required parameters for the corresponding equivalent 

s 

orthotropic plate. Thus 
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D = 


a(N x +1 ) 


^x^x 4 ®yl'y 4 
— n — m + — v n 
pV >s4 


and 


la 


a(N x +1 ) 


S” A 


p p 


[> + m y] 


(3.19) 


( 3 . 20 ) 


( 4 ) Square orthotropic grid having beams equally 
spaced in x and y directions ; 


In this special case of an orthotropic grid, the 
sides a and b of the grid are equal and the expressions 
given by (3.19) and (3*20) reduce to 
N +1 


D = 




5 \ ^c x x 


S_.lv m 4- + 


E y I y n j 


and 


h = 


VI 

-P n a 


c 


m x + m y ") 


( 3.21 ) 


( 3 . 22 ) 


Thus the analogy formulation for an qrthotropic 
grillage is complete. To validate this, a number of cases 
are considered in the next chapter. 


B - Isotropic grid 

An isotropic grid is a special case of the ortho- 
tropic grid when the spacing, material and geometric pro- 
perties of the beams are same in both the directions. 
Correspondingly an isotropic plate is a special case of 
orthotropic plate in which material properties have no 
preferred direction. Thus for an isotropic grid 


^x “y “ “b’ ^x -1 ^ ’ E x ~ ®y ~ 


(3.23) 
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ct „ b t\ — T) tj "n __ P~ 

HyTT - I x+ 1 ’ X - y - - - i - 12(1 .^2) 

From (3*15) and (3*17), one thus gets the flexural rigidity 
and thickness h for the equivalent isotropic plate as 


D. = 

l 


EfoLfo (N x +1 ) ^mVa 4 + nVb 1 ^ 
b £ m^/a^ + n^/b^-') 2 


(3.2*) 


2 (N x +1 ) 


(3.25) 


The special case of an isotropic grid is that with square 
boundaries and the above equations reduce to 


D. = 

l 


(N x +1 ) (m^ +n^ ) 


a (m^ + n^)‘ 


( 3 . 26 ) 


2 (M x +1) m b 


(3-27) 


For a square isotropic grid non-dimensional parameters 
flexural rigidity p and thickness parameter , q are defined 


¥b N x 


±I£_ 

2m b w x 


(3.28 


(3.29) 


Substituting expressions for and h from eq. (3*28) 
and (3.29) in eqs. (3*26) and (3.27) one obtains 


k 1+ 

, . 1 x m-tn^ 

P “ N x ( m 2«2)2 

* = (1 + 5T ) 


(3-30) 


(3-31 ) 



Note that ncodimensional flexural rigidity parameter p 
depends upon the order of the mode of vibration whereas 
the non dimensional thickness parameter q is independent 
of it . 


Similar dimensionless parameters, tension para- 
meter p^ = Sa/N^T and thickness parameter q^ = h a 2/ 
2A c aN x were defined for cable network"* , where 
A c = Area of cross -section of each cable. 

S = Tension per unit length of the equivalent membrane 
T = Tension in each cable. 

A 

These parameters were found to be equal and given by 

P F = q f , = 1 + 1/N X (3-31 'b' 

3 .1 *3 Grillages clamped on the boundaries : 

An orthotropic plate fixed on the edges is not 
exactly solvable but approximate methods have been used 
for its solution. Raleigh method has been used for rec- 
tangular orthotropic plate whereas results for square 
isotropic plate are known using Raleigh-? - * ^ method. For 
establishing analogy, expression for frequency derived by 
Raleigh's energy method is used here. 


The admissible function used for plate analysis 
is obtained from beam functions X^Cx), Y n (y) of the grill- 

27 

age and the displacement function w(x, y,t) is taken as 




B mn V x > V y) Sin V* 


(3-32) 


w (x,y,t) 
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The beam functions (eigen function) for c lamped - 
clamped beam are 


X r (x)=:Co S h if - Cos if - g r (Sinh f|^ - Sin ) (3.33) 

r = 1,2...... 

where e p are constants and their numerical values for 
each set of functions are given in Table-A. 


TABLE-A 


r 


e r 

i 

0.9825 

4.73 

2 

1 .0007 

7.8532 

3 

0.5359 

10.9956 

4 

1 .000 

14.1371 

5 

0.9999 

1 7.2787 

6 

1 .0000 

20.4203 

r > 6 

1 

(2r+1 )7? /2 


In Raleigh's energy method, one is mainly inte- 
rested in fundamental frequency and uses only X-j(;;) «vr\d 
Y-|(y). The admissible function for the plate, clamped on 
the edges, can thus be written, using eqs. (3*32) and 
(3-33), as 


r e m x 

B mn[ C ° Sh -S- 


! m x _ e m x 
a" - Cos "r 

im 5j[cosh ^ 


a * 

e n^ 

"ET " 


0 X 

l ( Sinh -SL 
m 


*E» 

‘H.l 


e n y 

Cos 


g. ,(Smh -v- 


¥)1 
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8 

Hearmon using Raleigh energy method, gives the expression 
for the fundamental frequency of the orthotropic plate, 
clamed on the edges, as 


P 


mn 


-/ 


i 


_fph 


A + n Al + 2H fafn (Oj, -2)(e n -2)l 
=* 7 h 1 * a? t? J 


a" 1 


• ... ( 3 - 35 

for, m = n = 1 . 

Let the flexural rigidity parameter for such 
a plate he defined as 

D oF - D x + 211 (e m" 2 ) ^ e n~ 2 ' 

'Then the frequency equation can be rewritten as 


(3*36) 


P mn = y D oF/fph m = n = 1 (3-37) 


A special case of the c lamped -c la mo oyrtLo^o^qi'c 
plate is an isotropic plate, clamed on the edges, for which 
Y^ung 2 ? found frequencies in different modes using Raleigh 
Ritz method. The calculation were carried out for a 36 
term series based on taking both m and n equal to 1,2,3 ,4, 

5 and 6. Young's re.sults for a square isotropic plate are 
reproduced in Table B. 

TABLE-B 

Mode NcuT T 2~~ 3 ~ 6 g 



35-99 73-41 108.27 131 .64 132,25 165-' 
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The rectangular isotropic plate b?s been solved 
2 12 14 - 17 

by various authers* *' ’ and their results (upto fifth 

frequency) are shown in the figure 3*3* The value of the 
parameter P mn b 2 y /^/D i , D ± being the flexural rigici^- 

ty parameter for isotropic plate can be read directly for 
a known a/b ratio. 


A Orthotropic grids: 


Proceeding in the same way as for simply supported 
case, the energy expressions for the grid and plate are 
found. Strain energy tr and kinetic energy X of the grid are 
derived and their expressions are 


T~ N 

~ ¥ B mn ^mn^^^x^x ( e m/ a ^ G m H n 

+ Eyly (e n A)) 3 G n hJU 7 J (3'3B1 

and 


Iv g 4 B mn '^mn 

Cos 2 

P t f 
unn | 

' m x a P 

'e~ G P H n 


mb 


N_ •* 



+ _I 

G1 n 

h/J 

- 

(3.3?) 

e n 




Q 

Whereas strain energy 0 Up and kinetic energy Xp of the 
plate are 


U 


- B 2 ^ Sin" 
P mn 2 


P t 
mn 


e 4 


e 4 
+ D — 2_ 
7 IT? 


2 H e m e n ( e m 2) (e n -2) 

P?2 


(3.40 


where, H = +2 
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K — -jL M p 2 B 2 Cos^ "d f a b 

P P, n .mn mn UU£> u „ o 


P 8 P *mn mn 

. N x n 

the constants Hn* , IL„ y 


) t ~ h_ Qi n-l . ( o 'Ll'] 

mn e e m n ’ u *^ 1 ' 

m n 


L n > 5 G m » G1 m are given below. 



C ° sh N,,+1 3 ~ Oos tJ 7~7T 3 ‘ ( Slr * tr.*i 


Ny 3* 

H /‘ = X 


X [cosh 

i=1 


N v+ 1 


- Sln ^ ) ] 2 


W x +1 " C ° S N x +1 g n (Sinh 
- Sin J±1 )1 2 

N x +1 J 


V 1 


( 3 . 42 ) 


(3.43) 


^r ^+S r ) (0.5 Sinh 2 e p + 2 Cosh e Sin b_) 


r\ 

+ ( 1 - S r ) (0. 5 Sin 2 e r + 2 Sinh e r Cos e y )+2 e, 
g r (Cosh 2 e p + 4 Sin e p Sinh e r - C: 


• (3- 


G1 p = (1+g r ) (0.5 Sinh 2 e r - 2 Cosh e r Sin e r ) 


2 

+ 0 ~g r ) ( 0 .5 Sin 2 e r - 2 Sinh e r Cos e r ) + 20 ^ 

- g r (Cosh 2 e r - if Sin e Sinh e„ - Cos 2 eJ 


(3M) 


r takes all the integervalues of m and n. 


2 


The constants G r ^ G1 p are calculated for ous values ol 
r and the results are given in Table 1 . It is to be not—" 
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that these do not depend upon the order of the grid that is 

N x 

number of beams N x or Ny, in the grid. Constants H n and 
N 

H m J are calculated for various order grids and depend upon 
N x or Their values are given in Table-2. Replacing 

expression in the brackets of equation (3 3+0 ) by D 0 jp 


u p = t ab d oF b2 Sln2 Pmn* 


(3*6) 


Equating maximum strain and kinetic energy expressions of 
the grid and the plate one gets the flexural rigidity pa- 
rameter and the thickness h, of the equivalent ortho- 
tropic plate and these are 


E X I X (a®) 3 G m H n x + Zyly ^) 3 G n ^(3*7) 


h - 2 e m e n 

* -fp ab 01 m 01 n 


m a N„ mb 

— G1 H n + -X. G1„ H 
p m a P n 

. hn u n 


'm J 


(3*8) 


To know the equivalent plate completely, one needs expr-^ 
ssions for D^, D and H. But here - one expression oq. 
(3*36) is available which only correlates these constants 
therefore it is not possible to determine the values of 
these constants individually. However, as for as the pro- 
blem of determining fundamental frequency is concerned 
one needs expression for D 0 jp as a whole. Thus it is suffi- 
cient to calculate value of D 0 p only. A special case of a 
clamped rectangular orthotropic grid is a clamped square 
orthotropic grid for which the required parameters are 
given as : 
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D oF = ^[Vx e m 3 G m V + Vy ®n 3 G n 


and 
h = 


2 e m e n 


Jp a G1 m G1 n 
B Isotropic Grids; 


r 


N 

% G ^m H n 


-m 


x 


N 

m y G ^n ^m 

en 


7 
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An isotropic grid is a special case of an ortho- 
tropic grid and for this the flexural rigidity parameter 
Dj_ and thickness h are directly obtained from expression 1 
(3.36), (3 • 1 +7) and (3.48) 


D i 


and 


Vb 

2 ab 


[ (jb 3 G m a 


r e m^~ e n + 0 emQn ( e m ~ 2 ) ( e n “2 ) *1 

L a 1 * T? a^b 2 J 


(3.5D 


h = 


Snb 

Jp 


I- 


N 

e n H n 


x 


N 


b G1 


+ 


e m H m 


y 


n 


3* G1 


(3-52) 


m 


In case of a clamped square isotropic grid, these para- 
meters take the form : 


N x N y 

n _ Vb f e m 3 G m H n + e n 3 G n H m y 

G i ~ 2 a 


[ 


e m 5- + e n 4 + 2 e m e n ( e m -2 ) ( e n -2 ) 


and 
h = 


^b 


a 


N x N 


G1 


n 


G1 


(3.53 


(3.5*+) 


m 


Thus the desired unknown parameters of the equivalent 
plate are known. The frequency calculated from these, giv 
the frequency of the grid. 


•* 
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3.2 SUMMARY OF AgSUhTS; 


The analogy proposed states that the frequency for 
the grillage is obtained by substitution of flexural rigidi- 
ty parameter and thickness of the equivalent plate in the 
frequency equation of the analogous plate. The results 
obtained for simply supported and clamped grillages are 
summarised below: 

A - Frequency equation for grids S.S. on the edges; 


As per the proposition simply supported grid is 
to be replaced by an equivalent plate with same boundary 
conditions. The proposition leads to frequency equation 
due to the fact that exact solution of plate with simply 
supported edges is known. For grid with different spacing 
and material properties in x and y directions, frequenev 
equation is 


Pmn ~TT 2 


E x Ix Ll EAT ), 

T~ ( N x + ^ ) m + (N v +1 ) n 

a3 b3 J 

m x a (N x +1 ) + rriyb (Ny+1 ) 


(3.5-53 


For orthotropic grids having uniform beam spacing in both 
x and y-direction, above equation reduces to : 


P. 


mn 


TT 




E V I 


'x x x m 


4 


3T 


+ h 

"b^ 


h 


(m 


■x + V 


( 3 . 56 ) 
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For a rectangular isotropic grid the frequency equation 



For square isotropic grid, S-S. on the edges, 


frequency equation reduces to 



(3 . 58 -a) 

(3. 58-Id) 


B Grillages clamped on the edge 


As mentioned earlier, no closed form solution is 
known for the frequencies of plates clamoed on the edges. 
However equation (3*37) can be used to give an approximate 
result for the fundamental frequency of an orthotropic 


plate. Thus the fundamental frequency of the orthotropic 
grids is obtained by substitution of D qF and h (equations 
3*^-7 and 3*^8) in equation 3*37* 



(3-59) 



3 *+ 


or 



For a square orthotropic grid, the fundamental frequency 
reduced to 


G1 1 G 1 e 1 2 


N* N v 

^X^X^I + Fyly H"| 

~ 5 * N~“~~ 

m X H 1 +m y H 1 


( 3 . 60 ) 


2 1 2 14- 1 7 27 

For isotropic plates several authors’ ' ’ ’ ’ „„ ...have, 

calculated the value of 'the parameter f defined as 


-• •• - __ •• mn /”D • ( 3 - 61 ) 

^ J~ v h iF 

m,n = 1,2 . , , 

for various modes and for various side ratios a/b of the 


plate. The flexural rigidity parameter and the thickness 
h, of the equivalent plate, derived in equations 3.5"1 and 
3 . 52 , are substituted in eq. 3*61 so as to give the fre- 
quencies of the grid, clamped on the edges.. 


The frequencies of the grillages are obtained by 
proposing the plate analogy. This proposition needs va- 
lidation. The proposition will be verified by taking vari- 
ous case studies. The details of the procedure are given 


in the next chapter 
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CASE STUDIES 


In chapter III, frequency equations for grids with 
two different boundary conditions, simply supported and 
clamped, have been derived from the proposition of the plate 
analogy. To validate the proposition, a number of grillages, 
mentioned below, have been solved by the finite element 
method and the results are compared with that of the pro- 
posed equivalent plate. 

A Grids, simply supported : 

(i) Orthotropic grids 

(a) Rectangular 

(b) Square 

(ii) Isotropic grids 

(a) Rectangular 

(b ) Square 

B Grids clamped; 

(i) Orthotropic grids 


(a) Rectangular 

(b ) Square 
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(ii) Isotropic grids 

(a) Rectangular 

(b) Square 

General treatment for cases (A) and (B ii) is 
given whereas case (B i) is dealt for its fundamental fre- 
quency only. 

4.1 FINITE ELEMENT METHOD : 

The Finite Element method is a well taiown techni- 
que for structural analysis and is not discussed in details. 
However its relevence to the present problem is discussed 
in the following steps: 

4.1.1 Basic element -of the grijf . 

The beam element between the two joints of the grid 
is fatten as the basic element.. As mentioned earlier, only 
bending vibration of the grillage have been considered 
and therefore the basic element of the grid has four degrees 
of freedom shown in figure 4-1 . Out of these two degrees 
of freedom are associated with the transverse deflection 
(V^ and V^) whereas the remaining two are associated 
with the rotations ( and . 

4.1.2 Numbering system of the joints : 

To denote the degrees of freedom some numbers 
have been assigned to these.. The numbering system for 
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the grid, S.S. on the boundaries, is shown in figure 4- -2 
and for the grid, clamped on the edges, is shown in 
figure 4- -3* Each internal joint is identified by a dis- 
placement coordinate and two rotations in xz and yz planes. 
Thus each internal joint is represented by three numbers of 
which, 

(1 ) First indicates the degree of freedom associated with 
the transverse deflection* 

(2) Second indicates the degree of freedom associated 
with the slope of the beam in x-jj plane- 

(3) Third denotes the degree of free tjcci associated with 
the slope of the beam in y-z plane. 

In case of S.S. grids (Fig* 4- -2), the joints on the boun- 
daries have only one degree -of freedom associated with the 
rotation whereas the joints of the grids clamped on the 
boundaries have no degree of freedom. 


4-. 1.3 Stiffhess and Mass matrices of the grid : 

The mass and stiffness matrices M and W 

« 

respectively of the frid is obtained by synthesis of mass 
and stiffness of the constituting element as per standard 
procedure. Let the np.ss'and stiffness matrix of the ele- 
ment be H 
is written as 


and 


It 


then the force deflection relation 


14- 
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0*.l) 


whe^S-^ and are shear forces and bending moments at the 
joint i and 


M - 


Vb 


V 1 


v 2 

© 2 

12 

61 

-12 

61 

61 

i + l 2 

-^1 

- 21 2 

-12 

-61 

12 

-61 

61 

21 2 

-6 1 

t - l 2 


(b.2) 


Correspondingly one can write the mass matrix (with cubic 
displacement curved) for the element as 


V20 


Vi 

©1 

V 2 

&2 

156 

22 1 

5+ 

_1 3i 

22 

b i 2 

131 

-3 1 2 

5V 

131 

156 

-221 

-13i 

-3i 2 

" 22 1 

^l 2 


(^•3) 


1 being the length of the basic element. 

V ..1 .4- Determination of frequencies : 

The differential equation of free oscillation of 
the discretized system is 


[m] 


q V + 


t l = 0 

* 


( U - 3 +) 



Having obtained |^MJ and l K J , one solves the equation 
( 4 . 4 ) for its eigen values and corresponding eigenvectors 
by Matrix Iteration method^ . 

A compute# program, refer Appendix C, is given to 
determine the frequencies and mode shapes of the rectangu- 
lar orthotropic grid. The programm furnishes as many 
eigenvalues as desired but gives mode shapes upto second 
only. 

4.2 LIST OF EXAMPLES : 


The list of grids considered for the dynamic 
analysis is given in Table C - The frequencies of all these 
grids are given in Appendix -B. The mode shapes (upto 
2nd mode) have also been found for all these cases. 

Define non dimensional frequency parameter 


U = 



A b ^ b a 1 

Vb 


0.25 


Q+-5) 


where a^ is the spacing of the beams. It is calculated 
for 3 x 3 isotropic grid from the proposed analogy method 
and its results are compared with those obtain from earlier 
methods^ refer Table- 3 * 
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4 . 3 DISCUSSION OF THE RESULTS AND VALIDATION 
OFTHE PROPOSED ANALOGY : 

The results of frequencies of different grids 
considered in the present work, refer Table -C, are given in 
Appendix-B. For simply supported grids, Tables 6, 7 and 
9, it is clear that for a square grid of specified dimen- 
sions and beam properties, the natural frequencies converge 
to those of the equivalent plate as the order of the' grid 
increases. This implies that the equivalent plate for all 
the grids of different orders with same dimensions and 
beam properties is a unique one. 

It is also observed that in some cases the results 
obtained by Finite Element method are higher than by 
proposed analogy whereas in other cases these are lower 
(Refer Tables 4 to 9). One of the reason for the erratic 
performance of the Finite Element method in these examples 
is the fact that compatibility conditions are not satisfied 
at the joints (Assumption of torsionless grid), The lack 
of compatibility induces excessive flexibility, while the 
assumed deformation pattern represent constraints which 
lead to excessive stiffness of the individual elements* 

The combination of elements is either too flexible or too 
stiff depending upon which effect predominates. 

The results obtained for different grids, given 
in Appendix-B, clearly show that the frequencies obtained 
by the proposed analogy are comparable with those obtained 



^3 


by- Finite Element method. Also table-3 shows the closeness 

i 

of the non dimensional frequency parameter u found by the 

i . i > 

proposed analogy method to that obtained earlier by other 
authors3 ,i+, 6’ 7 } 23 # , 

i 

Timoshenko proposed intutively the analogy of an 
orthotropic plate for static analysis of the beam grillage. 
According to his analogy, mass and flexural rigidity of 
the beams are assumed to be evenly distributed over the 
pitch of the grid. The proposition thus can be extended 
for dynamic analysis. Frederick and Falgout^ calculated 
the non dimensional frequency parameter based upon ortho- 
tropic plate theory and showed closeness • of his results 
to those obtained by discrete analysis. The results 
obtained in the present work are equally close to those 
of Finite Element method and those obtained by earlier 
research workers. However unlike Timoshenko's proposition 
of analogy based on intution, the present analogy is based 
on the rationale of energy equivalence and holds good for 
grids with any number of beams. 


Based upon Timoshenko's proposition following 
relations for the analysis of a torsionless isotropic 
grid are obtained. 


Fulv) 

DjL = — (N_+1 ) 

cl 


and 


aT b A b n x 


(4.6) 


h 


0+.7) 
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Ren ten derived, equilibrium equation for static 
analysis of an orthotropic grid as 


^x^x "(frw , c ^ J ' x 


+ - y ) ~& 4 v 

a 1 Sx 2 5y 2 


+ 


E I 
~y y 

a„ 


^W 


- 9xy 


. 8 ) 


For a torsionless isotropic square grid, contribution due 
to torsional rigidities are neglected and the equation 
reduces to 


Vb 


(4.9) 


fc 1 nx- by 4 ' ^ 

Timoshenko's proposition gives the following form of equa- 
tion of motion for transverse vibration of the equivalen- 
ts otropic square plate 


D. (^r + --) - _ V y, ~b 2 w 


(4.10) 


The frequency equation of (4.10) for simply supported 
boundary condition is 


P 


mn 


= TT‘ 


D i 


"h a 


4 


(m 4 + 


n 4 ) 


(4.11 ) 


Vp 


Substituting and h from eq. (4.6) and (4.7) in eq. 
(4.11), one gets the frequency equation for a torsionless 
square isotropic grid, simply supported at the edges as 


mn 



Vb 


=7r J afb Ab a 1 


.1+ 


O+Jp) (m* + h ) (4.12) 

X 



Note that as N is very large the frequency equation 
(4-. 12), approaches to that obtained hy proposed analogy, 
refer eq. (3*58-b). 

In the present formulation of plate analogy it 
has been assumed that the equivalent plate behaves as a 
real plate whereas Timoshenko’s proposition violates the 
real model of the plate in bending, refer eq. (4.10). 

Note eq. (4.10) is not of the form of plate in transverse 


vibration. 
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CONCLUSION 

In the present work a plate analogy of the 
torsionless beam grillage for its dynamic analysis has 
been proposed. The analogy has been proposed to arrive 
at simple procedure for dynamic analysis of the grids 
because discrete approach, such as Finite Element method^ 
is time consuming. The analogy proposed in this work is 
presented on the rationale of energy equivalence whereas 
Timoshenko's orthotropic plate analogy is proposed on the 
intuticn that for very large order of the beam grillage it 
approaches an orthotropic plate. 

The proposition of the plate analogy has been 
validated by c cm paring the results of the frequencies 
obtained by the proposed analogy method with those found 
by Finite Element method. These results also compare well 
with those obtained by earlier research workers. Note that 
the analogy proposition is easily formulated for the grill- 
age with simply supported edges due to the fact that the 
equivalent orthotropic plate with S.S. edges is exactly 
solvable. For orthotropic beam grillage with fixed ends 
frequency equations for 1st mode is derived using admissi- 
ble function and the proposition is validated. However 
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for isotropic beam grillage, clamped on the edges, 
results are obtained for higher modes as well which compare 
reasonably with Finite Element method. 

SCOPE FOE FURTHER WORK: 


Since the analogy proposition simplifies the 
method for dynamic analysis, it is proposed that further 
work should be done to validate the analogy based on 
energy equivalence to the cases of grids ^ith torsion. 

It will be quite interesting tb find out how far a pro- 
position can be formulated for non orthogonal beam grillage 
which may be more of academic interest then of engineering 
applications. 
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APPMDIX-A 


A -a Energy expressions f or ort h ot r o pic orthog ona 1 
g rid work. S.S. on al l edges : ' 


Taking the modal functions of beam to be obtained 
frcrn that of plate at known x and y, as mentioned in chap- 
ter 3, the mode shape for ith beam lying in x direction, 
can be written as (Refer eq. 3*2)* 

Vi = B Sin ISH2 Sin HUi. Sin p t (4-1) 

A mn a N +1 mn 

X 

The strain energy of a beam in bending is written as 


U- 


b 


Vb 


Xo 


( 


?> 2 w >2 

3x 2 


dx 


(A-2) 


Therefore the strain energy of the ith beam lying in x 
direction is 


U 


Exlx 


xi 



From equations (A-1 ) and (A-3) 


(A-3) 



_ E X I X 
2 


C f B ( f ^) 2 

j ^ mn a 


Sin 


rtfri 

N x+ 1 


Sin Sin 
a 



Vi 


^mn 


(SS) 


Sin 


2 mi 
Nt+1' 


Sin c 


P t 
Mnn 


du 


Sirr S2S dx 


3 


\ 


\ * - 




or 


52 


U 


i = (a?) 4 sii, 2 mi Rin 2 


»2 ^ Sin- Pmn t [a/2] ( AJf) 

The strata energy u xb of all beams lying along x-axis is 
obtained by summing the ab< 
to N. 


'ove expression over i, from i = i 


x 


U. 


xb ~ 


-^x^x 

2 


=4 (F) 4 [a/2] Sin 2 


N. 


x 


p mn t Y Sin* 1 _J2Er 

Si %+i 


or 


E.I.. 

U xb = -* X B 2 ^ 
xv 2 mn l a ; 


(a/2) [_£|fjsin 2 p m t (A.-?) 


Similarly one can write the strain 
lying along y-axis. Thus 


en 


er gy Uy_ b of* all beams 


^vl-i 


U , = Z p2 /mrr^ ,, , .ri r +n o 
yk 2 mn ^ET - ) (b/2) — -J Sin 2 p^ n t 


(A-6 ) 


Therefore the total strain 
grillage i 9 


energy U^, of the orthotropic 


or 


U g U xb + U yb 


(A-7) 


U - - mn 71^ ci- 2 r 4 

w ~ — — - J l ’ Qiv\ w VN J_ l 1 -‘-'-y- JL -v- y 

m — g - J L - (N 


- g - -g=-- • ■ Sin- p ffin t 


The kinetic energy express! 


a3 '"x +1) * 4 ^(V)3 

* * • (A —8 ) 

;ion for the transverse vibration 

of a beam is given as 11 
1 

^ = i j” m b t 2 dx where w = dw/dt 


(A-9) 
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Therefore the kinetic 
x -direction is 


energy 0 f ith beam lying in 



(A-10) 


Substituting w 

a 


from (A-1 ) in above expression, 


one gets 


8x1 = i S p 4 *4 Sin 2 SEE Sin 2 

O a Nx+1 

Since m x , the maarf/unit length of the beam, 
therefore, 


Cos 2 p t dx 
■unn UJS - 

( A— 1 1 ) 
is constant, 


Sin M2 


%i = l 2 PL B| n Sin 2 Mi cos 2 p t f 

■«x+» mn J 

*** (A-1 2) 

Further simplification gives:- 

^ = V 4 >4 S “ 2 8&C0.2 P mn t {a/2] (A-1 3) 

Thus the kinetic energy K* of all beMS) ^ 

14 2 o 

xi = 2“ %n B mn [a/2] Cos 2 p^t Sin 2 ^ 

i=1 N x+ 1 


N. 


K 


x 


"xb = 


i=1 


or 

^h = f p m n b 4 L a/2 J c °s 2 P mn t [£|fj (a-14) 

Similarly, one can write the kinetic energy K h of all 
beams lying m x -direction. Thus 7 


^=f4B^Ch/aJ cos 2 p mn t 


[ 


V 


(A-1 5) 


Therefore the total kinetic 
orthotropic grillage is obtained by 
1 3) and (a-14) : - 


5 ^ 

energy K of the 

o 

adding expression 


s 


(A- 


K g ~ K bx + ^ 




(A-16) 


Hence 


K ~ — rp f>2 p 2 , 

g 8 Pmn mn Cos P t 

u mu mn 


(V < N X +1 } + m y b (N y +1 >2 


A b ^erg_2Lexpress icns f 


(A-1 7) 




Strain energy expression for the bending of 
orthotropic plate is given by 2 * 1 -; 


a b 


3V 


Up= MH Dx ( ^ )2 + 2d i & & 


+ D_ (5^)2 

dy £ 


^ ^ ~7 


+ %.& 2 ' 


] 


dxdy 


CA-18) 


Substituting w, from ( A 

U 


P 


o o 


^ ) in (A-1 8) one gets" 

~ Sin 2 SEr S in 2 t 
b mn 


a b 

= * i J [a B m n <jpr Sin 2 as Sir, 2 


+ D :y B nn ( f ^ Sin 2 SS Sln 2 ngy Sln 2 ^ 

+ 21 Bmn m 2 n 2 Sin 2 p t c 2 mnx 2 
a^b^ mn s Cos^ 


a Cos^ I2S 
a b 


■J dxd 7 


where , H = d 1 + 20 


(A-1 9) 


xy 


* * * 



or 
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u. 


P = i Sin 2 P m t r m h % (ab) + Be 


+ 2 m2 n2 H 




Since ab Ap _ Area of the plate 
Hence 


,7 r ‘ f 

Sln2 W ’ 

r > 
[ D x ~ 

+ 2H 

si 7 

aV J 

L a 


+ D £L 

y ^ 


The kinetic energy expression for the transverse 
Of the plate is given below 11 


K 


a b 

j 5 M p^ dxd y 


o o 

Since M p , the mass per taait 
Therefore: 


area of the plate is 


M t 


a b 


5 I Vmn s ±n 2 Sin 2 S^Cos 2 

n r\ >> ^ 


K P = ^ 

p 2 J ^nn Jam 
o o 


or 




M^ab 


8 


^mn *mn 2os ^ P„..t 


mn 


(ab) 

(A-20) 


(A -21 ) 
vibration 

( A-22) 

constant 

P mn t dxdy 
(A -23 ) 


(A-2+) 
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TABLES OF THE HESUHPS 
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TAB LB-1 

Values of the constants G r and G1 p for the 


grids 

clamped on the 

edges. 

r 

G r 

Gl r 

i 

9 .46 

9-46 

2 

15.62 

15.62 

3 

29.87 

29.87 

4- 

27.27 

27.27 



Values 

TABLE-2 

W Y 

of the constant H r for 

vari ous 



values of 

r and Ny 


r 

V 2 

! 2 $ 

li 

U> 

II 

525 

»y = 5 

1 

3.053 

4.012 

5.004 

6.002 

2 

4.004 

4.136 

4.997 

5.950 

3 

0.3511 

6.076 

6.014 

7.119 


1 .982 

0.618 

7.221 

6.259 



TABLE -3 

Non dimensional . . , „ „ 

frequency parameter u.for a 3 x 3 square isotropic grid, S.S. on the edges 
m = 30 x 1 0 6 psi, = * 284 -, =3 in.^, 1 -^ = 0.25 in.^", a = 50 in. 
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T ABLE 4 -a 

Frequencies (cps) of a 2x3 orthotropic rectangular grid S.S. 
at the edges. 

E x = Ey = 30x10°psi , ^> x =^y = .284, A x = 3 in^, Ay = 1 .875 in^, 


B x = 

= 0.25 

in*, I = 0.088 in 1 * 

,a =50 incfy , a/b 

= 1 .332 

Mode 

No. 

Ellin gt on° 

Propo- 

Finite Element sed 
method method 

B 1 1 


53-7 

41 .78 

41 .7 

B 21 


84.8 

119.02 

117.8 

B 1 2 


153.0 

123.31 

124 .4 

B 22 


198.2 

167.85 

167.0 

B 31 


209-0 

258.62 

260.0 



TABLE 4-b 

t 


Frequencies 

(cps) 

of 3 x 4 orthotropic rectangular 

grid, 



S.S. at the edges. 





a/b = 1 . 25 



Mode 

No. 

6 

Ellingt on 

Finite 

Element 

method 

Proposed 

method 

B 1 1 


47.5 

39.19 

39.1 

B 21 


81 .25 

110.15 

110.1 

B 1 2 


154.2 

118.13 

117.5 

B 22 


178.9 

i 57 .il 

156.0 

B 31 


190.02 

235.43 

241 .5 
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TABLE 5 -a 

Frequencies (cps) of a2x3 orthotropic rectangular grid, S.S. 

at the edges . 

E x = 30x1 0 6 psi, Ey = 12 x 10 6 psi, ^ x = .23+, ^> y = .1, 

A x = Ay = 3 in. 2 , I x = Iy = 0.25 in. 4 a = 50 inch, a/b =1.332 


Mode No. 

6 

Ellingt cn 

Finite 

Element 

method 

Proposed 

method 

B 11 

56.3 

4 7 .' 4-3 

*+7.3 

b 21 

88.9 

131 .24 

130.15 

B 1 2 

168.1 

14-2.76 

142.6 

B 22 

210.5 

139.53 

139.2 

B 31 

230.5 

286.87 

28+. 2 


TABLE g-b 


Frequencies (cps) of a 3 x 4- orthotropic rectangular grid, 

S.S. at the edges, 
a/b = 1 .25 


Mode No. 


Ellington 6 Finite Proposed 
Element method 
method 


B 11 

52.5 

44.31 

44.2 

B 21 

86.2 

123.07 

128.1 

B 1 2 

153.5 

129.97 

129.5 

B 22 

192.3 

177.44 

176.8 

B 31 

225.5 

286.18 

281 .0 
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TABLE 6 -a 

Frequencies (cps) of 2x2 orthotropic square grid, S.S. 

at the edges. 

E x = e 7 = 3«x10 6 psi, =y y = 0.284, A x = 3 in 2 , 

Ay = 1 .875 in 2 , I x = 0.25 in 4 ", I = 0.088 in 1 *', a = 50 inch 


Mode No. 

Ellington^* 

Finite 

Element 

method 

Proposed 

method 

B 11 

33.5 

33.^5 

33^0 

B 1 2 

72.2 

73.50 

73.90 

B 21 

112.8 

116.04 

116.00 

B 22 

129.4 

135.01 

133.20 


TABLE 6-fc 

Frequencies (cps) of 3 x 3 orthotropic square grid, 

S.S. at the edges. 


Mode No. 

Ellingt on.6 

Finite 

Element 

method 

Proposed 

method 

B n 

33.0 

3 3* 1 43 

33.^0 

B 12 

73.6 

73.92 

73.90 

B 21 

115.2 

114.65 

116.00 

CM 

OJ 

163*0 

134.17 

133.20 

b 13 

147.2 

1 52 .41 

155.00 
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TABLE 6-c 

Frequencies (cps) of 4 x 4 orthotropic squire grid, 

S.S. at the edges. 


Mode 

Mo. 

Ellington^ 1 

Finite 

Element 

method 

Proposed 

method 

B 1 1 


33-25 

33-42 

33-40 

B 1 2 


73-7 

73.97 

73.90 

B 21 


115.00 

116.28 

116.00 

CM 

C\J 

pq 


133.00 

133.91 

133.20 

b 13 


155.00 

155. 16 

155.00 



TABLE 7- 

:§ 


Frequencies 

(cps) 

of a 2 x 2 square orthotropic grid, 



S.S. on the 

edges . 


= 30x10 6 psi, E 

y = I2xl0^psi, 

x = .284, 

'h 

11 

* 

S* 

A X = ^ = 3 

in B , 

I X = Iy = 0.25 

in* - , a = 

50 inch. 

Mode 

No, 

Ellingt on^ 

Finit e 

Element 

method 

Proposed 

method 

B 11 


37.00 

37.32 

37.25 

B 1 2 


84.2 

84.58 

85-7 

B 21 


125.0 

123.13 

127.2 

B 22 


144.5 

150.90 

148.2 
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TABLE 7-6 

Frequencies (cps) of a 3 x 3 square ort ^otropic grid, 

S.S. on the edges. 


Mode 

No. 

Ellin gt on^ 

Finite 

Element 

method 

Proposed 

method 

B 1 1 


37-2 

37.30 

37.25 

B 1 2 


85-0 

35.51 

85-7 

B 21 


126.5 

127.31 

127.2 

B 22 


14-7.2 

149.74 

143.2 

B 13 


172.5 

1 75-64 

182.0 



TABLE 7-c 



Frequencies 

(cps) 

of a 4 x 4 square orthotropic grid, 



S.S. on the 

edges . 


Mode 

No. 

Ellingt orfi 

Finite 

Element 

method 

Proposed 

method 

B n 


37.23 

37.29 

37.25 

B 1 2 


85.52 

35.66 

85.7 

B 21 


127.0 

127.73 

127.2 

B 32 


147.9 

149 .40 

143.2 

B 13 


178.2 

1 30.40 

132.0 
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TABLE 8 -a 

Frequencies (cps) of a 2 x 3 isotropic rectangular grid, 

S.S. on the edges. 

\ = 30x10 6 psi A b = 3 in 2 1 ^ = 0.25 in 4- , = .284, a = 50in. 

a/b = 1 .332 


Mode 

No. 

Ellingt on^ 

Finite 

Element 

method 

Proposed 

method 

B 1 1 


52.6 

52*87 

52.8 

B 21 


112.5 

112.47 

113*2 

B 1 2 


179*9 

1 78.36 

1 36.2 

B 22 


204.9 

224.51 

211 .5 

B 31 


223*2 

- 

237*0 



TABLE 8- 

b 


Frequencies 

(cps) 

of a 3 x 4 isotropic rectangular grid. 



S.S. on 

the edges 




a = 5 ° in 

a/b = 1 . 25 


Mode 

No. 

Ellingt on^ 

Finite 

Element 

method 

Proposed 

method 

B 1 1 


47*7 

43.08 

43.05 

B 21 


110.9 

112.16 

• 111 .3 

B 1 2 


163*0 

- 

164.2 

B 22 


189*5 

- 

192.5 

B 31 


22 a .5 

- 

237*0 
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TABLE 9 - a 

Frequencies (cps) of a 2 x 2 isotropic square grid, 

S.S. at the edges. 

E fe = 30x1 0 6 psi, A b = 3 in 2 , I b = 0.25 irr, = 0.284, 

a = 50 inch. 


Mode No. 

Ellingt on^ 

Finite 

Element 

method 

Proposed 

method 

B 11 

36.6 

36.67 

36.64 

B 1 2 

103-6 

106.34 

106.8 

B 21 

IO 3.6 

106.34 

106.8 

B 22 

142.0 

143.31 

146 .6 


Frequencies 

(cps) 

TABLE 9-b 

of a 3 x 3 isotropic square 

S.S. at the edges. 

grid, 




Finite 

Proposed 

Mode 

No. 

Ellin gt on° 

Element 

method 

method 

B 11 


36.59 

36.65 

36.64 

B 1 2 


106.00 

1 06 . 79 

106.8 

b 21 


106.00 

106.87 

106.8 

B 22 


145.10 

147.18 

146.6 

B 13 


241 .00 

231 .08 

234.6 
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TABLE 9-c 

Frequencies (cps) of a 4x4 isotropic square grid, 

S.S. at the edges. 


Mode No. 

Ellingt on^ 

Finite 

Element 

method 

Proposed 

method 

B 1 1 

36.59 

36 .65 

36.64 

B 1 2 

106 .4 

.1 06 . 71 

106.8 

B 21 

106 .4 

106.83 

106.8 

B 22 

146.2 

144.67 

146.6 

B 13 

232.2 

233.96 

234.6 


TABLE 10 -a 

Fundamental frequencies (cps) of rectangular orthotropic 

grids clamped on the edges. 

E x = E y = 30 x 10 6 psi, A x = 3 in 2 , Ay = 1 .875 in 2 , 

I x = 0.25 in 4 , I = 0.088 in 4 , V x = = 0.284, 

a = 50 in 


Order of the 

grid 

a/b ratio 

Finite Element 
method 

Proposed 

method 

2 x 

3 

1 .332 

94.76 

93.4 

3 x 

4 

1 .25 

88.89 

88.2 

4 x 

5 

1 .2 

85-72 

85-5 
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TABLE 10-b 

Fundamental frequencies (cps) of square orthotropic grid 

clamped on the edges 


Order of 
grid 

the Finite Element 
method 

Proposed 

method 

2x2 

76.06 

75*5 

3x3 

75*86 

75*5 

4x4 

75*80 

75*5 

5x5 

75*73 

75-5 


TA BLE 11 -a 

Frequencies (cps) of a 2 x 3 isotropic rectangular grid 

clamped on the edges 

E^ = 30 x 10 6 psi , A b = 3 in 2 , 1^ = 0.25irf, ^ ='284, 

a = 50 in, a/b = 1 .332 


Mode No. 

Finite Element 
method 

Proposed 

method 

B 11 

1 1 9 * 79 

11 7-2 

b 21 

193-05 

189*0 

B 1 2 

270.02 .. . 

259*0 

B 22 

309.02 

303*0 
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TABLE 1 1 -b 

Frequencies (cps) of a 3x4 isotropic rectangular grid 

clamped on the edges 
a/b = 1 .25 
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table i 2 -a 

Frequencies (cps) of a 2 x 2 square isotropic grid 

clamped on the edges 

Eb = 30 x 10 6 psi, A b = 3 in 2 , 1^ = 0.25 in 4 , ^ = .28+, 

a = 50 in. 


Mode No. 

Finite 

method 

Element Proposed 

method 

B 11 

83-41 

82.8 

B 1 2 

162.04 

160.8 

B 21 

1 62.04 

160.8 

B 22 

233.55 

225.2 


table 12 -b 

Frequencies (cps) 

of a 3 x 3 

square isotropic grid 


clamped on 

the edges 

Mode No. 

Finite 

method 

Element Proposed 

method 

B 11 

83.18 

82.8 

B 1 2 

1 70.89 

1 70.0 

B 21 

1 70.89 

1 70-0 

B 22 

'231. #0 

225.2 

B 13 

292.10 

291 .1 
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TABLE 12 -c 


Frequencies (cps) of a 4 x 4 square isotropic grid 

clamped on the edges. 


Mode No. 

Finite Element 
method 

Proposed 

method 

B 1 1 

83.11 

82.8 

B 1 2 

1 71 -98 

1 70.2 

B 21 

1 71 .93 

1 70.2 

B 22 

229.90 

225.-2 

B 13 

316.04 

302.0 


TABLE 12 -d 


Frequencies (cps) 

of a 5 x 7 square 

isotropic grid 


clamped on the edges 

Mode No. 

Finite Element 
method 

Proposed 

method 

B 11 

83-09 

82.8 

B 1 2 

1 72.18 

1 71 .1 

B 21 

1 72-18 

1 71 .1 

cm 

CM 

PQ 

229.43 

226.1 

b 13 

320.89 

319.0 
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TABLE 13 -a 

First mode shape of 2 x 2 orthotropic grid, S.S. on the 
edges, with same dimension as given in Table 6 -A 


Finite 1 .000 

13.801 
Element -0.500 
-1 .000 

Method 

1 .000 
0.500 
0.502 
-1 .000 

1 .004 
-0.502 

1 3 • 801 

1 3 .'801 
-1 .004 
-0.500 

0.500 

1 .004 
- 0.502 

0.502 
13.801 
-1 .0C4 

Proposed 1 .000 

1 .000 

1 .000 

1 3 • 783 

0.500 

■ 0.500 

13*783 

0.500 

-0.500 

-1 .000 

1 .000 

13-783 

Method -0.500 

0.500 

13.783 

-0.500 

-0.500 

-1 .000 

-1 .000 

-1 .000 







TABLE 13 -fa 



First mode shape 

of 3 x 3 orthrotropic 

grid, S. 

. S . on the 

edges 






Finit e 1 . 000 

1 .414 

1 .000 

1 .001 

11 .253 

0.707 

0.708 

15.921 

1 .000 

0.000 

1 1 .258 

0.707 

Element -0.708 

-1 .001 

1 .416 

15.921 

0.000 

1 .001 

22-516 

0.000 

0.000 

15.921 

0.000 

-1 .001 ' 

Method -1 .416 

1 .001 

11.258 

-0.707 

0.708 

15.921 

-1 .000 

0.000 

11 .258 

-0.707 

-0.708 

■41 . 001 

-1 .000 

-1 .414 

-1 .000 





Proposed 1.000 

1 .414 

1 .000 

1 .000 

11 .254 

0.707 

0.707 

15. 916 

1 .000 

0.000 

11 .254 

0.707 

Method -0.707 

-1 .000 

1 .41 4 

15.916 

0.000 

1 .000 

22.508 

0.000 

0.000 

15.916 

0.000 

-1 .000 

-1 .414 

1 .000 

11 .254 

-&.? 07 

0.707 

15. 91 6 

-1 .000 

0.000 

11 .254 

-5.707 

-0.707 

-1 .000 

-1 .000 

-1 .414 

-1 .000 
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TABLE 13-c 

Second mode shape of a 2 x 2 orthotropic grid, S.S. at the 

edges . 


Finite 1 .000 

-13.471 

Element -0.500 
-1 .000 

Method 

-1 .000 
-0.500 
-0.877 

1 .000 

1 *754 
-0.877 
-13 .4 72 

13-472 

1 -754 
0.500 

0.500 

1-754 

-0.877 

-0.877 

13.472 

1 .754 

Proposed 1 .000 

-1 .000 

2.000 

13.750 

0.500 

-0.882 

-13.750 

-0.500 

-0.882 

2.000 

2.000 

13.750 

Method -0.500 

-C.882 

-13.750 

0.500 

-0.882 

2.000 

-1 .000 

1 .000 






Second 

mode shape of a 

TABLE 13-d 

3x3 orthotropic 
edges . 

grid, S 

. S . on the 

Finite 

1 .000 

0.000 

-1 .000 

1 .269 

1 1 .1 75 

0.707 


0.000 

0.000 

0.000 

-1 .369 

-11 .1 75 

-0.707 

Element 

0.000 

1.396 

1 -936 

15.305 

0.000 

0.000 


0.000 

0.000 

-1 -936 

-15.805 

0.000 

0.000 

Method 

1 -936 

1 .369 

11 .175 

-0.707 

0.000 

0.000 


0.000 
-1 .000 

-1 .369 
0.000 

-11 .1 75 

1 .000 

0.707 

0.000 

1 .369 

Proposed 1 .000 

0.000 

-1 .000 

1 .423 

11 .231 

0.723 


0.000 

0.000 

0.000 

-1 .423 

-11 .231 

-0.723 

Method 

0.000 

1 .423 

2.000 

15.850 

0.000 

0.000 


0.000 

0.000 

-2.000 

-15.850 

0.000 

0.000 


2.000 

1 .423 

11 .231 

-0.723 

0.000 

0.000 


0.000 
-1 .000 

-1 .423 
0.000 

-11 .231 

1 .000 

0.723 

0.000 

1 .423 



APPENDIX -C 
C OMPOTER 


PROGRAMME 



oo ooooooooooo 




DYNAMIC .ANALYSIS OF ORIHCTROPIC 

RECTANGULAR GRID SIMPLY SUPPORTED ON THE BOUNDARIES. 
THIS PROGRAM GENERATES THE MASS AND STIFFNESS 
MATRICES OF THE GRID AND CALC JUTES EIGEN VALUES 
AND EIGENVECTORS OF THE GRID. 

S REPRESENTS STIFFNESS MATRIX. 

G REPRESENTS THE MASS MATRIX. 

AR=AREA OF C ROS SECT I ON . W=W EIGHT DENSITY. 
BR=BREDTH,H=HEIGHT OF THE CROSSSCTION. 

E=ELASTICITY CONSTANT, DR=M0MANT OF INERTIA, 
X=LENGT., R=MASS DENSITY. 

common /bi /s(5o,5o) 

C OMMON/B 2/G ( 5 0 3 0 ) 

COMMON/B3/ND 
common/b4/C(.5o,5o) 

DATA E,BR,H, XI, W/30.0E+06, 3 - , 1 .,50.0,0. 2$f/ 

DATA El ,BR1 ,H1 , Wl/12, 0E+06, 3 • ,1 • , • V 
IR=ER*H**3/1 2. 

DR1=BR1*B±**3/12. 

AR=ER*H 
ARl=3!l*Hl 
R=V386.4 
R1 =WT. '386 .4 

NY REPRESENTS THE NO. OF BEAMS IN X -DIRECTION 
AND NX REPRESENTS NO. OF BEAMS IN ¥- DIRECT ION 
D01 000NY=2 ,4,1 

45 FORMAT (1 OX *01=* El5*5 ,*C2=*E1 ,*C3=*, E15-5, 

*c4=*,bi 5 - 5 ) 

46 FORMAT (IX, *CRDER OF THE GRID IS=* ,I2,*BY* ,12) 

NX=NY+1 
MS=NX+1 

X=X1 /FLOAT (MS ) 

Cl=E |£ DR/X**3 
C 2=R* AR* X/4 20 . 

ERINT45,C1 , C 2 , C 3 , C4 
C NA IS THE ORDER OF THE MASS AND STIFFNESS MATRICES 

NA=3*NY** 2+7*NY+2 
D01 1=1 ,NA 
D01 J=1 ,NA 
G(I, J)=0. 

1 S(I,J)=0 
D 021=1 ,NY 

S(I,I )=C1*4 .*X**2 
G(I, J )=C2**4.*X**2 

2 CONTINUE ' - 

NC=NA-NY 

D 031=1 ,NY 
ND=NC+i 

S(ND,ND)=C1*4 ,*x**2 
G(ND,ND)=C2*4 ,*X**2 



3 CONTINUE 
' Dd+I=1 ,NY 
J=3*I-2 

ne=nx+j 

S(I,NF)=-C1-*6«*X 
S(I,NF+1 )=C1*2.*X**2 
G(I ? NE)dC2*13.*X 
G(I,NF+1 )=-C2*3»*X**2 
k- CONTINUE 
D05J=1 NY 
1=3* J -2 

NI =3* NY* *2+3* NY +1 +1 
NJ=(5*NI**2+1 7*NY)/2-1 4J 
S(NI,NJ')=C1*6 .*X 
S(NI+1,NJ)=C1*2.*X**2 
G(NI,NJ)=-C2*13.*X 
G(NI+1 ,NJ)=-C2*3.*X**2 
5 CONTINUE 
NLJf*NX+3 
NI=3*NY-2 
D07K=1,NI,3 

NM=NX+K 

NN=NL+K 

S(NM,NN)=-C1*1 2. 
S(NM,NN+1 )=C1*6 .*X 
S(NM+1 ,NN)=-C1*6.*X 
SCNM+1.NN+1 )=C1*2.*X**2 
G(NM,NN)=C2*5if . 

G(NM,NN+1 )=C2*13**X 
G(NM+1,NN)=C2*13.*X 
G(NM+1 ,NN+1 )=-C2*3 .*X**2 
7 CONTINUE 
NP=3*NY+2 
I L=N Y-1 
D 06 1=1 ,IL 
D06K=1 ,NI,3 
D06I=1 ,2 
D06J=1 .2 
IM=NX+I 
IN=NL+J 
IQ=IM+L*NP 
IR=IN+L*NP 
IM=IM+K-1 
IN=IN+K-1 
NQ=NQ+K-1 
NB=NR+K-1 

s(nq 3 nr)=s(im,in) 

G(NQ,NB)=G(lM, in ) 
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6 


CONTINUE 

S(NX,NX)=C3*If .*X**2 
SfNX>NX+1 )=-C3*6.*X 
S(MjNX+3 )=C3*2.*X**2 

G(NX,NX)=C l f*l+.*X**2 

G(NX,NX+1 )=C4*'!3.*x 

G (NX -NX 4-3 ) =-Cb* 3 . *x* * 2 

NS=3*NY-5 

D08K=1 ,NS,3 

NT=NX+K 

SfNT ,NT )=2?+ .*C1 +£+ ,*C3 
SfNT , NT +3 )=-C3*1 2 . 
S(NT,NT+5)=C3*6.*X 
S(NT+1,NT+1 )=sC1*8 »*X**2 
S(NT42,NT+2)=C3*8.*X**2 
SiNT +2>NT+3)=-C3*6.*X 
S(NT+2 ) NT+5’)=£!3*2.*X**2 

+3 )=&■’ >C1 +24 .*G3 
S (NT -fif ,NT )=C1 *8 .*X**2 
S(NT+5,NT+5)=C3*8.*X**2 
G(NT,NT )= 3 12.*C2+ 3 12.*C4 
G( NT ,NT + 3 )=£! ! +*5 i + • 

g(nt,nt+5)=-c4*i-3.*x 

G(NT+1 ,NT+1 )=C2*8.*X**2 
G ( NT+2, NT +2 ) =C4*8 . *X** 2 
G(NT+2jNT+3 )=C 1 +*1 3 .*X 
G ( NT +2 , NT +5" ) =■ -C4-* 3 . +X 5 ** 2 

J >JJJ+3 ,NT +3 ) =31 2 .*C2 +3 1 2 .*c4 
G(NT -*4 ,NT -4+ )=C2*8 ,*X**2 
G (NT+5 ,NT +5 ) =C4* 8 .*3?=* 2 
CONTINUE 
NU=4*NY+2 

S(NU,NU)=C 3 *4 ,*X**2 
S(NU-1 , NU)=C3* 2.*X**2 
S(NU-3,NU)=C 3 *6.*X 

G CTJU , NU ) =04*4 ,*X** 2 
G(NU-1 ,NU)=-C4*3 .*X**2 
G(NU-3,NU)=-Clf*13.* X 
D09K=1 ,NY 
D 091=1 ,NP 
D09J=1 ,NP 

nv=i+k*np+ny 

NW=J +K*NP+NY 

MA=NY+I 

MB=NY+J 


S(N7,NW)=S(MA,MB) 

G(NV,NW)=$(MA,MB) 
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100§ 


CONTINUE 
DOI 0K=2, NA 
MC=K-1 
D010J=1 M3 

S(K,J)=S(j k) 
C (KjJ )=G(Uj k) 

CONTINUE 
PRINTS, NY, NX 
CALL BINV 
CALL MAT ML 

call eignv 
CONTINUE 
STOP 
END 


£IBFTC SUE1 


STOROOTINE OT BIOT CAIC:JL ^ rE THE IIIVEilsE 0F AM AT BIX 
COMMON/BI /A( 50, 50) 

C0MM0N/B3/N 

DIMENSION INDEX (50 3 ) ’ 

S-f A V) iUmCS (IEOW ’ jicW A(ICOUJM,JCOLTM),(iM ; lX ) I, 
DO 20 J=1 N 
INDEX (I, 3 ^ 

DO 550 1=1 N 

SEABCH FOB PIVOT ELEMENT 

AMAX=o.n 

D01 05J=1 N 

1P5.60 

TV ^SC^OO 80 

“ AVS ( A( J > K) 5 ^85 , 1 00, 1 00 

IC OIIJM=K 
AMAX=ABS(A( J, K) ) 

c ontinue 

CONTINUE 

™5aS 3 o^ IM(IOOffiM >3 >+i 

INDEX (1,2) =IC DLUM 

S T SSf G H0WS T0 PDT PIT ® elsmskis 

DO 200 L=1 N ’ ’ ° 

S WA P=A ( IBOW . L ) 

A(IR0W L)=AOICOLUM L) 

a(icohjm,l)=swap ’ 

divide pivot bow by pivot element 
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31 0 PIVOT =A( ICQLCJM, ICQLUM ) 

330 A(IC0LUM, IC OLUM ) =1 .0 
340 D0350L=1,N 

350 A( IC OLUM, L ) =A( IC OIXJM, L ) /PIV CT 

C REDUCE NON -PIVOT ROWS 

380 D0550L1=1,N 

390 IF (LI -IC OLUM )400, 550,400 
400 t=a(li,icolum) 

420 A( LI , IC 0LUM)=0 . 0 

430 DC45*0L=1% N 

450 a(li ,l)=a(li,l)-a(icolum,l)*t 

55O CONTINUE 

C INTERCHANGING COLUMNS 

600 DO71 01=1 ,N 
610 LAN+1 -I 

620 IF( INDEX (L,1 ) -INDEX ( L, 2) )630, 71 0,630 

630 JROW=JNDEX( L, 1 ) ' ' 

640 JC OLUM=INDEX ( L, 2) 

650 D0705K=1,N 

660 SWAP=A(X. JROW) 

6 70 ‘ A(K,JHOW)rrMK/UQ-QUJM) 

700 A( K, JC 0LUMT=SWAP '' 

705" CONTINUE ^0 RETURN 

710 CONTINUE / SND 

THIS SUBROUTINE IS TO CALCULATE THE 
MULTIPLICATION OF TWO MATRICES 
SUBROUTINE MATML 
COMM ON /B1 /SC 50,50) 

C 0MM0N/B2/G ( 50, 50 ) 

C OMMON /B3/NA 

common/e4/d(5o,5o) 

D01I=1,NA 
D01 J=1 ,NA 
D(I,J)=0. 

D01 K=1 ,NA 

D(I,J)=D(I,J)4S(I,K)*G(X, J) 

1 CONTINUE 
RETURN 
END 

£ IB FTC EIGNV 

THIS SUBROUTINE CALCULATES EIGEN -VALUES AND 
EIGEN -VECTORS 
SUBROUT I 'IE EIGNV 
C0MM0N/B3/M 
C01-JMCN/E4/C (50,50) 

DIMENSION EIGN(9),Y(50),X(50),Y1 (5$ z(5o) 
DIMENSION X2(50),X1 (50),D (' 1 j5°) 

4 FORMAT (12, *TM. EINENVALUE=* , El 6 .8,1GX,*FBE- 
' UENC Y=* , FI 0 . 2,1 OX , 14 ) 

5 FORMAT ( 1 2, *TH. EIGENVECTOR* ) 

6 FORMAT (8F10 f 5) 
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EQUIVALENCE (Y1,Z) 

E=0 .00061 
D01 1 0IV =1 ,5 
JK=IV 
ITER=0 

IF(M.EQ.I) Y(1 )=1 .0 
IF(IV-M)116,117,116 

1 16 N=M-IV+1 
N1 =N-1 
D0100 1=1 ,N 

100 Y1 (I )=1 . 0 . 

99 KC=0 

C AND EIGENVECT OES ARE CARRIED CUT 

DO 101 1=1, N 
Y(I)=0.0 
DO 10T J=1,N 

101 Y(I )=Y(I )■+€ (iiU)* Y1 (J) 

R=Y(1 ) 

DO 9 © 1=1 ,N 

Y(I )=Yf I)/R 

IF(ABS(Y(I )/Y1 ( I ) — 1 .0 ) -E) 90,90,102 

102 KCil 
90 CONTINUE 

ITER=ITSR+1 

IF(ITER.GE.300) GO TO 103 
IF(KC) 103, 103, 1 0+ 

1 04- D 0115 1 = 1 , N 

119 Y1(I)=Y(I) 

GO TO 99 

C UPIO THE STATEMENT 113, THE SIZE OF THE C 

C MATRIX IS REDUCED BY ONE 

- 103 EIGN(IV )=R 

DO 109 1=1 ,N1 

109 X(I)=-Y(I+1) 

DO 121 J=2,N 

121 Z(J)~C(1,J) 

IFCIV.GT.1 ) GO TO 8 
D07I=1 N 

7 X2(I)=Y(I) 

D01 J=1 ,M 

1 D(1,J)=C(1 ,J) 

8 CONTINUE 

DO 112 1=2, N 
DO 112 J=2,N 

112 C(I -1 ,J-1 )=0(I, J)+Z(J)*X(I-1 ) . 

GO TO 118 

117 EIGN (XV )=C(1 ,1 ) 



118 FREQ =1 -0/(SQRI(EIGN(lV))*6.2832) 
PRINT V, IV, SIGN (IV ), FREQ, ITER 
IFCIV.EQ.1 ) GO TO f+1 
IF ( IV . GT .2) GO TO 110 
XI (1 )=o. 

IC =M-1 
D 021=1 - IC 

2 X1(I+1J=Y(l) 

SUH=0. 

D 091=1 ,M 

9 SUM=SUM+D(1 ,I)*X1 (I) 

F=(EIGN(2)-EIGN(1 ))/SUM 
D03I=1 ,M 

3 X(I)=X2(l)+F*X1 (I) 

19-1 PRINT 6, (Y(I),I=1 ,M) 

110 CONTINUE 

RETURN 

END 


GENTRY 



